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ABSTRACT 

We compare observed clustering of quasars and galaxies as a function of redshift, mass, luminosity, and 
color/morphology, to constrain models of quasar fueling and the co-evolution of spheroids and supermassive 
black holes (BHs). High redshift quasars are shown to be drawn from the progenitors of local early-type 
galaxies, with the characteristic quasar luminosity L* at a given redshift reflecting a characteristic mass of the 
"active" BH/host population at that epoch. In detail, the empirically calculated clustering of quasar "remnants" 
(knowing the observed clustering of the original quasars) as a function of stellar mass and/or luminosity is 
identical to that observed for early-type populations. However, at a given redshift, the active quasars clus- 
ter as an "intermediate" population, reflecting neither "typical" late nor early-type galaxies at that redshift. 
Comparing with the age of elliptical stellar populations as a function of mass reveals that this "intermediate" 
population represents those ellipticals undergoing (or terminating) their final significant star formation activ- 
ity at the given epoch. Assuming that quasar triggering is associated with the formation/termination epoch 
of ellipticals predicts quasar clustering at all observed redshifts without any model dependence or assump- 
tions about quasar light curves, lifetimes, or accretion rates. This is not true for spiral/disk populations or the 
quasar halos (by any definition of their ages); i.e. quasars do not generically trace star formation or halo assem- 
bly/growth processes. Interestingly, however, quasar clustering at all redshifts is consistent with a constant host 
halo mass ~4x 10 12 /z _1 Mq, similar to the local "group scale." The observations support scenarios in which 
major mergers dominate the bright, high-redshift quasar population. We demonstrate that future observation of 
quasar clustering as a function of luminosity can be used to constrain different fueling mechanisms which may 
dominate AGN populations at lower luminosity and/or redshift. We also show that clustering measurements at 
Z = 3 - 6 will be sensitive to the efficiency of feedback or "quenching" at these redshifts. 

Subject headings: quasars: general — galaxies: active — galaxies: evolution — cosmology: theory 



1. INTRODUCTION 

In recent years, it has become clear that essentially 
all galaxies harbor superma ssive black holes (BHs) (e.g., 
iKormendy & Richstonell9 95). the masses of which are corre- 
lated with many properties of their host spheroids, including 
luminosity, mas s (Magorrian et al. 1998), velocit y dispersion 
dFerrarese & Merrittll2000t iGebhardt et al.ll2000l) . concentra- 
tion and Sersic index (iGraham et al.ll200lh . Further, compari- 
son of the re lic BH mass density and quasar luminosity func- 
tions (OLFs:ISoltanll982tlSalucci et al.ll999HYu & Tremainel 
2002tlMarconi et al.l2 004: Sha nkar et al. 2004; Hopkin s et al l 



2007bl th e cosmic X-ray backgroun d (Elvis et al.1 120021: 



Ueda et al. 2003; Cao 



2005; Hopkins et al.ll2007bl). and relic 



BH Eddington ratios (Hopkin s. Hernquist. & Naravari 2005) 
demonstrates that the growth of BHs is dominated by a short- 
lived phase of high accretion rate, bright quasar activity. This, 
together with the similarity between the cosmic star formation 
history (e.g., Hopkins & Beacom 2006) and quasar accr etion 
history (e.g.. iMerloni et all 120041; iHopkins et al.ll2007bl) . re- 
veals that quasar activity (with associated BH growth) and 
galaxy formation are linked. 
A number of theoretical models have been proposed to 
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explain the evolution of these popul ations with redshift, and 
their co rrelations with on e anoth er jKauffmann & H aehneltl 
[20051 ISomerville et alJ l200Tt IWvithe & Loebl 120031; 
Grana tcTet al.1 12004 IScannapieco & Ohl 120041; iBaugh et al.l 
20051: iMonaco & Fontanotl 120051; ICroton et al I | 2006f 
Hopkins et al.l l2006bllcl l2007at iBoweretaTT bood: 
iDe Lucia et al.l 120061: iMalbon et ail 120061). In many of 
these models, the merger hypothesis dToomrdl 1977b provides 
a physical mechanism linking galaxy star formation, mor- 
phology, and black hole evolution. According to this picture, 
gas-rich galaxy me rgers channel large amounts o f gas to 
galaxy centers (e.g jBarnes & Hernquist ! 1 99 U[l99 6). fueling 
powerful starbursts (e.g., [M ihos & H ernquisn 1 19961) and 
buried BH growth (Sanders et al. 1988) until the BH grows 
large enough that feedback fro m accretion rapidly unbinds 
and heats the surrounding gas (Silk&Rees 1998), briefly 
revealing a bright, optical quasar ( Hopk ins et al 1 12005 al) . 
As gas densities and corresponding accretion rates rapidly 
decline, an inactive "dead" BH is left in an elliptical galaxy 
satisfying observed correlations between BH and spheroid 
mass. Major mergers rapidly and efficiently exhaust the 
cold gas reservoirs of the progenitor systems, allowing the 
remnant to rapidly redden with a low specific star formation 
rate, with the process potentially accelerate d by the expulsion 
of gas by the quasar ( Springel et al. 2005al). 

Recent hydrodynamical simulations, incorporating star for- 
mation, supernova fee dback, and BH growth and feedback 
dSpringel et al.l 12005b) make it possible to study these pro- 
cesses dynamically and have lent support to this general sce- 
nario. Mergers with BH feedback yield remnants resem- 
bling observed ellipticals in their correlations with BH proper- 
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ties (|Di Matteo et al.l2005l). scaling relations (Rob ertson et alj 
2006b), colors (Spring e! et al.ll20 05a). and morphological and 
kinematic properties dCox et alj|2006l) . The quasar activity 
excited through such mergers can account for the QLF and 
a wide range of quasar prope rties at a number of frequencies 
(Hopkins et al. 2005a, 2006b), and with such a detailed model 
to "map" between merger, quasar, and remnant galaxy popu- 
lations it is possible to show that the buildup and statistics of 
the quasar, red galaxy, and merger mass and luminosity func- 
tions are consistent and can be use d to predict one another 
dHopkins et al.ll2006d.l2007all2006et) . 

However, it is by no means clear whether this is, in fact, 
the dominant mechanism for the triggering of quasars and 
buildup of early-type populations. For example, the associa- 
tion between BH and stellar mass discussed above leads some 
models to tie quasa r activity directly to star formation (e.g., 
iGranato eT al. 2004), implying it will evolve in a manner trac- 
ing star-forming galaxies, with this evolution and the corre- 
sponding downsizing effect roughly independent of mergers 
and morphological galaxy segregation at redshifts z < 2. Oth- 
ers invoke post-starburst AGN feedback to suppress star for- 
mation on long timescales and at relative ly low accretion r ates 
through, e.g., "radio-mode" feedback (Croton et al. 2006|). In 
this specific case, the "radio-mode" is associated with low- 
luminosity activity after a quasar phase builds a massive BH 
(i.e. qua sar "relics"), but it is possible to con struct scenar- 
ios (e.g., ICiotti & Ostrikerll200U iBinneyl 120041) in which the 
same task is accomplished by cyclic, potentially "quasar- 
like" (i.e. high Eddington ratio) bursts of activity, or in which 
the "radio-mode" might be directly associated with an opti- 
cally luminous "quasar mode," either of which would imply 
quasars should trace the established "old" red galaxy popu- 
lation at each redshift. In several model s, a distinction be- 
tween "hot" and "cold" accretion modes ([Birnboim & Dekell 
120031: iKeres et al.ll2005t Ibekel & Birnboimll2006l) . in which 
new gas cannot cool into a galactic disk above a critical dark 
matter halo mass, determines the formation of red galaxy pop- 
ulations, essentially independent of quasar triggering (e.g., 
ICattaneo et~ail 2006, but see also Binney 2004). 

At low luminosities, (M B > -23, important at z < 
0.5), models predict that stochastic, high-Eddington ra- 
tio "Seyfert" activity triggered in gas-rich, disk-dominated 
systems will contribute increasingly to th e AGN luminos- 
ity function dHopkins & Hernquisj 120061) . with enhance- 
ments to these fueling mechanisms from bar instabilities 
and galaxy harassment. Indeed, the morphological makeup 
of low-luminosity, low-redshift Seyferts appears to support 
this, with increasing dominance of unperturbed disks at 
low luminosities seen locally (e.g., [KauffmannetalJ|2003b; 
iDong & De Robertisl 120061) and at low redshift z ~ - 1 
(ISanchez et al.ll2004t iPierce et all 12006). At high luminosi- 
ties (even at these redshifts), however, the quasar popula- 
tions are increasingly dominated by ellipticals and merger 
remnants, particularly those with y oung stellar populations 
suggesting recent sta r burst activity (iKauffmann et al.li2 003b ; 
Sanch ez et ail l2004t IVanden Berk et alj 12006; Be st et al.1 
2005i: lDong & De R obertis 2006), and even clear merger rem- 
nants ( San chez et al.l 120041) . Still, some models extend the 
observed fueling in disk systems to high redshift quasars, in- 
voking disk instabilities in very gas- rich high redshift disks 
as a p rimary triggering mechanism (IKauffmann & Haehneltl 
l2OO0l) . 

Clearly, observations which can break the degeneracies be- 
tween these quasar fueling models are of great interest. Un- 



fortunately, comparison of the quasar and galaxy or host lu- 
minosity functions, while important, suffers from a number of 
degeneracies and can be "tuned" in most semi-analytic mod- 
els. Direct observations of host morphologies, while an ideal 
tool for this study, are difficult at high redshift and highly in- 
complete, especially for bright quasars which dominate their 
host galaxy light in all observed wavebands. However, the 
clustering of these populations may represent a robust test 
of their potential correlations, which does not depend sensi- 
tively on sample selection. Critically, considering the cluster- 
ing of quasars and their potential hosts is not highly model- 
dependent in the way of, e.g., mapping between their luminos- 
ity functions or modeling their triggering rates in an a priori 
fashion. 

In recent years, wide-field surveys such as the Two Degree 



Field (2dF) QSO Redshift Survey (2QZ; iBovle et al. 
and the Sloan Digital Sky Survey (SPSS: lYorketal 



2000) 



2000) 



have enabled tight measurements of quasar clustering to red- 
shifts z ~ 3, and a detailed breakdown of galaxy clustering as 
a func tion of galaxy mass , luminosity, color, and morpholog y 
(e.g., iNorberg et aLll2002h IZehavi et al.ll2002t iLi et al.ll2006l) . 
These observations allow us to consider the possible trigger- 
ing mechanisms of quasars in a robust, empirical manner, 
and answer several key questions. Which local populations 
have the appropriate clustering to be the descendants of high- 
redshift quasars? How is the quasar epoch of these popula- 
tions related to galaxy formation? And, to the extent that 
quasars are associated with spheroid formation, are bright 
quasar populations dominated by quasars triggered in forma- 
tion "events"? 

In this paper, we investigate the link between quasar activ- 
ity and galaxy formation by comparing the observed cluster- 
ing of quasar and galaxy populations as a function of mass, 
luminosity, color, and redshift. In § 12 we compare the clus- 
tering of quasars and local galaxies to determine which galaxy 
populations "descend" from high-redshift quasar progenitors. 
In §|3] we consider the clustering of quasars as a function of 
luminosity and redshift, checking the robustness of our re- 
sults and presenting tests for the dominance of different AGN 
fueling mechanisms at low luminosities. § [4] compares the 
clustering of quasars as a function of redshift with that of dif- 
ferent galaxy populations at the same redshift, ruling out sev- 
eral classes of fueling models. §[5] further considers the age 
as a function of stellar and BH mass of these galaxy popula- 
tions, and uses this to predict quasar clustering as a function 
of redshift for different host populations. In §[6] we use these 
comparisons to predict quasar clustering at high redshifts, pre- 
senting observational tests to determine the efficiency of high- 
redshift quasar feedback. Finally, in §|7]we discuss our results 
and conclusions, and their implications for various models of 
quasar triggering and BH-spheroid co-evolution. 

Throughout, we adopt a WMAP1 (Dm, D a, h, erg, n s ) = 
(0.27 0.73,0.71,0.84,0.96) cosmology dSpergel et alj 
120031) . and normalize all observations and models shown to 
this cosmology. Although the exact choice of cosmology 
may systematically shift the inferred bias and halo masses 
(primarily scaling with a&), our comparisons (i.e. relative 
biases) are unchanged, and repeating our calculations for 
a "concordance" (0.3,0.7,0.7,0.9,1.0) cosmology or the 
WMAP3 results of Spergel et al. (2006) has little effect on 
our conclusions. All magnitudes are in the Vega system. 

2. USING CLUSTERING TO DETERMINE THE PARENT 
POPULATION OF QUASARS/ELLIPTICALS 
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At a given redshift Zu quasars are being triggered in some 
"parent" halo population. These halos, and by consequence 
the quasars they host, cluster with some bias/amplitude b(zd- 
The halos will subsequently evolve via gravitational cluster- 
ing, which in linear theory predicts their subsequent clustering 
at any later redshift zj will be given by 



b(z f )=l + 



D(zd 
D( Zf ) 



b{zd-\ 



(1) 



jFrvl fl99l IMo & Whitd fl996t iCroom et all l200lh . where 
D(z) is the linear growth factor. Thus, at z = 0, the halos 
which hosted the quasars at Zi will have a bias of b(0) = 
l+D(zi)[b(zd-ll. 

The quasar luminosity function at a given redshift has a 
characteristic luminosity L*. Given that quasars (at least 
those with L > L„) are typically obser ved to have high Ed- 
dington ratios A = L/L P m 0.3-0.5 dHeckman et alj|2004t 



Vester gaardll2004t iMcLure & Dunlopl l2004t iKollmeier et alj 
20051) . this L* reflects the characteristic mass of "active" 
BHs at that redshift, M BH « 3.0 A -1 x 10 8 M Q (L*/T0 13 L Q ). 
Direct observations of quasar Eddington ratio s/BH masses 
dMcLure & Dunlopll2"004t IKollmeier et alj|2005l) : limits from 
the X-ray backgroun d dElvis et al.| 120021: ICaol 12005)) and 
BH mass functions dSoltanl ll982UYu& 
Marconi et all l2004t IShankar et all 12004 



2007bl). rad io luminosity 



1J jZ-UU-t 

i 982t lYu"& Tremaind 1200 2 : 

Hopkins < 
functions ( Merloni 2004 



: et al. 
2004; 



McLure et a l. 2004), and "relic" low-lumin osity Eddington 
ratios (Hopkins, Hernquist, & Naravan 2005) all rule out the 
possibility that these BHs subsequently gain significant mass 
(> 10-20% growth) after their brief "active" phase, so the 
Mbh above is equivalent (within its errors) to the z = BH 
mass of these objects. Since the relationship between BH 
mass and host luminosity or mass is well-determined at z = 
(M BH = M^gai with fi » 0.001; Marconi & Hunt 2003, Hiiring 
& Rix 2004), knowing M Bli {z = 0) of a population implies, 
with little uncertainty, its z = host mass M ga i or luminosity. 

In FigureQ] we consider the clustering of quasars as a func- 
tion of redshift, evolved to z = 0. At each redshift where the 
quasar clustering b(z) is measured, we also know the char- 
acteristic luminosity L*. Here we adopt the bolometric 
determ ined in the observational compilation of Hopkins et al. 
d2007bl) . 

(2) 
(3) 



log(L* /L Q ) = k + ki £ + k 2 £ 2 + h £ 3 , 

e^iog(^) b ref =2], 



(k = 13.036; k\ = 0.632; k 2 = -11.76; k 3 = -14.25), but 
for our purposes this is iden t ical to adopting the fl-ban d 
or X-ray L* f r om lUeda et alj d2003l); ICroom etalJ d2004l) : 
Hasin ger et alj d2005l) ; [Richards et alj d2006a) and converting 
it to a bolomet r ic L» with a typical bolometric c orrection from 
IMarconi et all d2004l) : iRichards etafl d2006bl) : iHopkins et ail 
d2007bl) . Given the conversions above, we consider the im- 
plied characteristic BH mass and, assuming little subsequent 
BH growth, the corresponding z = stellar mass or luminos- 
ity in a given band ( here from the Mbh _ ^host relations of 
IMarconi & Huntll2003l) . Knowing how the bias of these halos 
evolves to z = (Equation [TJ, we plot the bias as a function 
of stellar mass, at z = 0, of the evolved quasar "parent" pop- 
ulation. We compare this with observed bias as a function of 
stellar mass or luminosity for both early and late-type galax- 
ies. Note that the lowest and highest-redshift bin s (z ~ 0.5 and 
Z ~ 2.5, respectively) in the IM vers et al.l d2007al) quasar clus- 
tering measurements are significantly affected by catastrophic 



redshift errors (owing to their considering all photometrically 
classified quasars, as opposed to just spectroscopically con- 
firmed quasars); we follow their suggestion and decrease (in- 
crease) the clustering amplitude in the lowest (highest) red- 
shift bin by ~ 20%; excluding these points entirely, however, 
has no effect on our conclusions. 

For reference, we show the characteristic Mbh correspond- 
ing to L* in the QLF (as defined above) as a function of 
redshift in Figure [2] Whether we adopt a direct conversion 
from the observed L* (Equation (01 ) with observed Eddington 
ratios, as above, or invoke any of several empirical models 
for the QLF Eddington ratio distribution, we obtain a similar 
Mbh- The figure illustrates the inherent factor < 2 systematic 
uncertainty in the appropriate Eddington ratios and bolomet- 
ric corrections used in our conversions. These uncertainties, 
however, are at most comparable to the uncertainties in quasar 
clustering measurements. Because of this, our conclusions 
and comparisons are not sensitive to the exact method we use 
to estimate the Mbh corresponding to L* at a given redshift. 

The comparison in Figure [T] is possible at any redshift, not 
simply z = 0. We repeat our methodology above at several 
Zobs = 0-1, evolving the bias to b{z bs)- The agreement with 
red galaxy clustering observed as a function of mass at each 
Zobs is good, at all z < 1 . At higher redshifts, small fields in 
galaxy surveys limit one's ability to measure clustering as a 
bivariate function of luminosity and color/morphology at the 
highest luminosities, where the relics of z ~ 2 - 3 quasars are 
expected. 

3. CLUSTERING AS FUNCTION OF LUMINOSITY AND DIFFERENT 
AGN FUELING MECHANISMS 

The comparison in Figure Q] has one important caveat. We 
assumed that measurements of quasar clustering at a given 
redshift are representative of a "characteristic" active mass 
Mbh oc of the QLF. In other words, quasar clustering 
should be a weak function of the exact quasar luminosity, at 
least near L*. If this were not true, our comparison would 
break down on two levels. First, it would be sensitive to 
the exact luminosity distribution of observed quasars. Sec- 
ond, if quasars of slightly different luminosities at the same 
redshift represented different BH/host masses (consequently 
making quasar clustering a strong function of quasar lumi- 
nosity), there would be no well-defined "characteristic" active 
mass at that re dshift. 

Fortunately, ILidz et all d2006l) considered this question in 
detail, and demonstrated that realistic quasar light cu rve and 
lifetime models like those of IHopkins et alj (2006b) indeed 
predict a relatively flat quasar bias as a function of luminosity, 
in contrast to more naive models which assume a one-to-one 
correlation between observed quasar luminosity, BH mass, 
and host stellar/halo mass. This appea rs to be increasingly 
confirm e d by direct ob s ervati o ns, with [A delberger & Steidel 
(2005a); Croometa l] d2005l); iMvers et alj d 2006l l2007al): 
iPorciani & Norberdd2006l) : ICoil et all d2006bl) finding no ev- 
idence for a significant dependence of quasar clustering on 
luminosity. 

Figure [3] explicitly considers the dependence of bias on lu- 
minosity and its possible effects on our conclusions. We plot, 
at each of several redshifts, the observed bias of quasars as 
a function of luminosity. For the sake of direct compari- 
son, all observations are converted to a bolometric luminosit y 
with the bolometric corrections from Hopki nset alJ d2007bh . 
The QLF break lumin osity at each redshift, estimated in 
Hopkins et al. (2007b), is also shown. The first thing to note 
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Fig. 1 . — Evolved clustering of quasar "descendents" (colored points) as a function of mass or luminosity, compared to clustering of early type (solid black 
points) and late type (open black points) galaxies of the same mass/luminosity. The measured clustering of quasars at a given z. (samples as labeled) is ev olved in 
linear theory to the given observed redshift, and plotted as of function of the relic host gala xy mass/luminosity . Galaxy clustering is shown at z = from | Li et al.l 
12006, color and mo rphologically-selec ted early types as squares and ci rcles, resp ectively), Zehavi et al. 12005, color-selected; stars), Percival et al. (2006, SDSS 
LRGs; triangles), Norberg et al. 12002, diamonds) ; at z = 0. 1 —0.4 fr om Shepherd et al. (2001, stars) and Brown et al. (2003, and references therein; diamond); at 
z = 0.4-0.8 from Meneux et al. 12006, square), Phleps et al. 12005, circle), and Brown et al. (2003, and references therein; diamonds); and at z = 0.8-1.5 from 
Meneux et al. ( 2006, square), and Coil et al. 1 2004, circles). Fitted lines show the best-fit bias of early type galaxies at each z as a function of mass/luminosity 
INorberg et al.l200ll;lTegmark et aU2004l ; ILi et alJ2006D . 

per logarithmic interval in relic mass. The integrand here 
defines the relative contribution to a given observed lu- 
minosity interval from each interval in Mbh- Given the 
BH-host mass relation, we can convert this to the rela- 
tive contribution from hosts of different masses M ga i (i.e. 
d0g(L)/dlogM ga i). In detail, we assume P(Mbh M ga i) is dis- 
tributed as a lognormal about the mean corre lation, with a dis- 
persion of 0.3 dex taken from observations dMarconi & Hunt] 
l200l iHaring & Rixl l200l: iNovak et al.ll2006h and hvdrodv- 
namic al simulations (Di Matt eo et al.l 12005: Rob ertson et al.l 
l2006al) . Calculating the bias for a given relic Mbh or M ga i and 
observed redshift as in §|2l we can integrate over these contri- 
butions to determine the appropriately weighted mean bias as 
a function of observed quasar luminosity, 

(b)(L, z) = f b(M s . dh z) ^ Q{ ^ dlogMgaj. (5) 

<j> Q (L) J dlogM ga i 

Although binning by both luminosity and redshift greatly 
reduces the size of observed samples and increases their er- 
ro rs, the observatio ns in Figure [3] confirm the predictions 
of iLidz et al.l d2006l) to the extent that they currently probe. 
To contrast, we construct an alternative "straw-man" model. 
Specifically, we compare with the naive expectation, if all 



is that the quasar observations with which we compare gen- 
erally sample the QLF very near L», so regardless of the 
dependence of bias on luminosity, our conclusions are not 
changed. We have, for example, recalculated the results of 
§ [2] assuming that the characteristic mass of active BHs is 
given by M Edd ((L obs )), where (L bs) is the mean (or median) 
observed quasar luminosity in each clustering sample in Fig- 
ure Q] and find it makes no difference (changing the compar- 
isons by <C Ict). 

We compare the observations with various theoretica l mod- 
els in Figure |3] The models of Hop kins et all d2005bl) define 
the conditional quasar lifetime; i.e. time a quasar with a given 
final (relic) BH mass (or equivalently, peak quasar luminos- 
ity) spends over its lifetime in various luminosity intervals, 
fg(L|MBH)- Since this is much less than the Hubble time at 
all redshifts of interest, the observed QLF (</>g[L]) is given by 
the convolution of tQ(L\M^n) with the rate at which quasars 
of a given relic mass Mbh are "triggered" or "turned on," 



<p Q (L)= / f e (L|M BH )<?XMBH)dlogMBH, 



(4) 



where 0(Mbh) is the rate of triggering, i.e. number of 
quasars formed or triggered per unit time per unit volume 
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■ Kollmeier et al. 2005 (Measured: BH Virial Relation) 

• Hopkins et al. 2006 (Simulations+Observed QLF) 

* Yu S Tremaine 2002 (Continuity+BHMF) 
Merloni et al. 2004 (QLF FP) 
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FIG. 2. — Characteristic "active" mass of BHs at a given redshift, i.e. the BH 
mass corresponding to (dominant at) L„ in the QLF. Black squ ares adopt the 
virial relation BH mass determinations of Kollmeier et al. 1 2005). Red circles 
fit the observe d QLF at each redshif t (Hopkins et al. 2007b) to quasar light 
curve models I Hopkins et al. 2006b). Blu e sta rs adopt the s implified con- 
tinuity model from Yu & Tremaine (2002) and Marconi et al. 1 2004) given 
the lUeda et al] 120031) QLF Blue dashed line I Merloni 2004) and derives the 
active mass distribution from the radio and X-ray "black hole fundamental 
plane." Black solid line is a fitted relation (Equations {5) & |SJ). The "active" 
BH mass is well-defined. Adopting different estimat ors does not significantly 
alter our conclusions or comparisons in Figuresm& ll II 

quasars were at the same Eddington ratio (so-called "light- 
bulb" models), i.e. if there was a one-to-one correlation be- 
tween observed luminosity and BH mass (and correspond- 
ingly, host mass), which produces a much steeper trend of 
bias as a function of luminosity and is significantly disfavored 
(> 4.5 cr; although from any individual sample the signifi- 
cance is only ~ 2 a). We also compare the predicted cluster- 
in g as a function of luminosity fro m t he semi-analytic model s 
of lKauffmann & Haehneltl (120021) and lWyithe & Loebl (120031) . 
which adopt idealized, strongly peaked/decaying exponential 
quasar light curves (i.e. Eddington-limited growth to a peak 
luminosity, then subsequent L oc exp(-f/fg)) and therefore 
yield similar predictions to the constant Eddington ratio "light 
bulb" model (and are likewise disfavored at > 4a). 

Figure |4]highlights the dependence of bias on luminosity in 
the observations and models by plotting the relative bias b/b* 
(where = b(L*)) near the QLF L„, more clearly demon- 
strating the observational indication of a weak dependence. 
Alternatively, we can fit each observed sample binned by lu- 
minosity at a given redshift to a "slope," 



b d(b/b*) 
b* dlogL 



(6) 



the results of which are shown in Figure [5] as a function of 
redshift, compared to the slope (evaluated at L„) predicted by 
the various models. At all redshifts, the observations are con- 
sistent with no dependence of clustering on luminosity, and 
strongly disfavor the "light-bulb" class of models (again, at 
~ 4a at z ~ 1-5 — 2). This confirms the conclusions of these 
studie s individually, partic ularly the most recent observations 
from iMvers eta l. (2007a ]) and the largest luminosity base- 
line observations from Adelberger & Steidel (2005a). The 
weak d e pendence predicte d by the models of Hopki ns et al.1 
(2006bT); lLidz et all {2006 ) provides a considerably improved 
fit, although even it may be marginally too steep relative to 
the observations. 

Galaxy clustering (and therefore, presumably, host halo 
mass) appear to be much more strongly correlated with galaxy 



luminosity or stellar mass (Figure [TJ than with quasar lu- 
minosity (at a given redshift); i.e. the weak dependence of 
bias on quasar luminosity appears to be driven by variation 
in Eddington ratios at a characteristic "active" mass. This is 
also supported by comparison of quasar lumi nosity functions 
and n umber counts, in a semi-analytic context (IVolonteri et al.l 
2006). We note that this is completely consistent with obser- 
vations that find similar high Eddington ratios for all bright 
quasars, as these are c onfined to L > (and indeed the 
iHopkins etail d2006bl) : iLidz et alj d2006l) model predictions 
do, at these highest luminosities, reproduce this and imply a 
steep trend of bias with luminosity). However, the relatively 
weak trend in clustering near and below makes our conclu- 
sions throughout considerably more robust, so long as the ob- 
served quasar sample resolves (true for all plotted points). 

Despite the detail of the models involved, the predictions in 
Figure [3] are all simplified in that they model only one mech- 
anism for quasar fueling. However, Hop kins & Hernquisl 
(2006) (among others) predict that at low luminosities, con- 
tributions from smaller BHs in non-merging disk bulges, 
triggered by disk and bar instabilities, stochastic accretion, 
harassment, or other perturbations, are expected to domi- 
nate the "Seyfert" population. We therefore repeat our cal- 
culation, but allow different fueling mechanisms in differ- 
ent hosts to contribute to different quasar luminosities , ac- 
cording to the m odels of Hopkins & Hernquist (200g) and 
ILidz et al] (120061). Because the "Se yferts" in this particular 
model iHopk ins & Hernquistl 120061) are generally less mas- 
sive systems at high Eddington ratio in blue, star-forming 
galaxies, they are less biased than merger remnants of sim- 
ilar observed luminosity. 

The inclusion of these populations in Figure [3] does not 
change our conclusions near L ~ L*. However it does in- 
troduce a feature, generally a sharp decrease in observed 
bias, at the luminosity where these secular fueling mecha- 
nisms begin to dominate the AGN population. This luminos- 
ity is typically quite low, L ~ 10 11 - 10 12 L© (corresponding 
roughly to luminosities below the classical Seyfert-quasar di- 
vision of Mb = -23). The only redshift at which the cluster- 
ing of such very low-luminosity AGN has been measured is 
z < 0.2, by which point massive, gas-rich mergers are suf- 
ficiently rare that the predi cted "Seyfert" population from 
Hopkins & Hernquistj (120061) dominates the merger-triggered 
quasar population at all luminosities, erasing the feature in- 
dicative of a change in the characteristic host population. 
However, it is possible that deeper clustering observations at 
Z ~ 1 -2 will eventually probe these luminosities, and test this 
prediction. 

Realistically, the luminosities of interest are sufficiently low 
that X-ray surveys present the most viable current probe, but 
with the small < 1 deg 2 field sizes typical of most surveys, the 
quasar autocorrelation function cannot be constrained to the 
necessary a ccuracy to distinguish the mod els. However, as 
proposed in Kauff mann & Haehneltl d2002l) . the AGN-galaxy 
cross-correlation presents a possible solution. For example, 
there are sufficient galaxies at z ~ 1 in the fields of surveys like 
e.g. DEEP2 or COMBO-17 (with field sizes ~ 3.5 deg 2 and 
~ 0.8 deg 2 , respectively) that the accuracy of cross-correlation 
measurements is limited by the number of AGN; considering 
the hard X-ray selected AGN samples in the CDF-N or CDF- 
S (withheld sizes -0.01-0.5 deg 2 ) fromz - 0.8- 1.6 would 
represent a fact or — 2 - 3 increas e in the number density of 
AGN over the ICoil et al.l d2006bl) sample in Figure [3] while 
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FIG. 3. — Quasar bias as a function of quasar luminos ity at each of severa l redshifts. Points show quasar observations from Figure[T](same style). Open black 
points add the local (z < 0.3) observations of IConstantin & Vogelevl <2006l SDSS LLAGN (LINERS +Seyferts); squares). IWake et all <200H SDSS; star), and 
IGrazian et al. 1 2004, AERQS; diamond), and the z ~ 1 cross-correlation measurements of ICoil et alj (2006b, SDSS+DEEP2; circles). Blue dotted line shows 
Z.*(z). These are compared t o vario us models described in the text (curves, as labeled). Models in black adopt the feedback-regulated quasar light curve/lifetime 
models from Hopkins et al. (2006a), others consider more simplified "light bulb" model light curves. 
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FIG. 4. — As Figure [3] but showing the relative quasar bias b/b„ 
(b* = b(L*)) as a fu nction of luminosit y near L* for the models and 
most well-constrained (Myers et al. 2007a) and largest luminosity baseline 
l Adelberger & Steidel 2005a) observations, to highlight the luminosity de- 
pendence and differences between the models. 



FIG. 5. — The best-fit dependence of quasar bias on luminosity 
(d(fo/fc*)/dlogL) from the observations in Figures[3]&[4](points), compared 
to the dependence expected from the models (calculated at L„). The observa- 
tions favor little or no dependence of clustering on luminosity. 

extending the AGN luminosities to a depth of ~ 3 x 10 9 L Q (~ 
10 42 ergs" 1 cm" 2 ). The measurement of the cross-correlation 
between observed galaxies in these fields and deep X-ray se- 
lected faint AGN at z > 0.5 does, therefore, present a realis- 
tic means to test the differences in these models at low lumi- 
nosities. The observation of a feature as shown in Figure [3] 
should correspond to a characteristic transition in the quasar 
host/fueling populations. 
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4. CLUSTERING OF DIFFERENT POPULATIONS 

In Figure[6]we compare the observed quasar bias and corre- 
lation length as a function of redshift with the expected clus- 
tering of quasar hosts, i.e. evolving the observed bias of BH 
(quasar "relic") hosts up from z = 0. A z = elliptical galaxy 
or spheroid of stellar mass M ga i has a bias b(M, z = Q) shown 
in Figure [1] and a BH of mass Mbh ~ 0.001 M ga \. For con- 
venien ce, we adopt the analytic fit to Z?(M ga i, z = 0) in lLi et all 
(HH), 

2>(M gaI ,z= 0)/K = 0.896 + 0.097 (M gaI /M„), (7) 

where here M * = 1.02 x 10 n M fn is the Schechter function 
"bre ak" mass (|Bell et al.l 120031) and b * = 1.2 for red galax- 
ies (IZehavi et alJ 12005b iLi et all 120061) . The progenitors of 
these systems therefore represent the "characteristic" active 
systems when the QLF characteristic luminosity L*(z) is given 
by L*(z) « ALEdd(0 001M ga i), i.e. when these BHs domi- 
nated the ~ L* quasar population and assembled most of their 
mass. Evolving the local observed bias of the z = spheroids, 
b(M g a\, z = 0), to this redshift with Equation (fl} yields the ex- 
pected bias that the quasar hosts (and therefore quasars them- 
selves) at this redshift should have, £>g(z). For future compar- 
ison, this is approximately given by 



b Q {z)~ l+0.0UD(zT l W 



5.70a-2.30je -3.35 r 1 



(8) 



where x = log(l +z) and D(z) is again the linear growth fac- 
tor. This expectation is plotted, with the ~ ler combined un- 
certainties from errors in the measured QLF L» and local bias 
b(M ga \, z = 0), comparable to the inherent factor < 2 system- 
atic uncertainty in the appropriate Eddington ratios and bolo- 
metric corrections. We also plot the corresponding correlation 
length ro; because measurements of this quantity are covariant 
with the fitted correlation function slope 7, we renormalize 
the models and observations to 



l Mpc( 



ro.fit 



8/r'Mpc 



7/1.8 



(9) 



This is similar to the non-linear <r^ L parameter (standard de- 
viation of galaxy count fluctuations in a sphere of radius 
ShT^Mpc, i.e. a$ measured for an evolved density field, see 
PeebleJ |1980b . and effectively compares the amplitude of 
clustering at 8/r'Mpc with a fiducial model with 7 = 1.8, 
minimizing the co variance. 

The expectation agrees well with observed quasar cluster- 
ing as a function of redshift (x 2 /v = 29.6/32, with no free 
parameters). For comparison, we plot the expected cluster- 
ing of halos of a fixed mass, M ha io ~4x 10 11 - lO 13 /z _1 M , 
det ermined in th e conte xt of linear collapse the ory follow- 
ing iMo & White! dl996l) . modified according to Shet h et all 
(2001) and with the po wer spectrum calculated following 
lEisenstein &Hul {1999) for our adop ted cosmology. As 
noted in most pre vious studies (e.g., iPorciani et al.l l2004t 
ICroom et al.l [2005), a constant host halo mass of a few 
10 12 h~ l Mq provides a surprisingly good fit to the trend with 
redshift. This empirical fit, is comparable to our expecta- 
tion from elliptical populations (best-fit halo mass 3.86 x 
10 12 h~ l M Q with x 2 /V = 28.9/31; of course, the exact best-fit 
mass depends systematically on cosmology, but this conclu- 
sion is robust). There is at most a marginal trend for the halo 
mass to increase with redshift (forMh a i oc (1 +z) , the best-fit 
k = 0.41 ± 0.45; corresponding to a ~ 50% increase over the 
observed redshift range). 



Note that the measurements shown are not all statistically 
independent, and the significance of this comparison will di- 
minish if we consider any single quasar clustering measure- 
ment. Figure [^demonstrates the same comparison, highlight- 
ing individual quasar bias measurements separately. How- 
ever, the previous agreement and our co nclusions are similar 
in all c ases. As discussed by the authors. IPorciani & Norbergl 
(2006) find a somewhat higher c lustering amplitude in their 
highest-redshift (z *~ 2) bin than Cro om et al.1 d20 05) study- 
ing the same sample ( a nd high er than IM vers et al.l d2006l) and 
lAdelberger & Steide l (2005a) w ho study independent sam- 
ples), but the significance of the IPorciani & Norbergl (120061) 
result is < 2a. 

That a constant halo mass fits the data as well as observed 
suggests that there may be a physical driver or triggering 
mechanism associated with these halos. It is suggestive that 
this corresponds to the "group scale;" i.e. minimum halo mass 
of small galaxy groups, in which galaxy-galaxy mergers are 
expected to proceed most efficiently. However, the redshift 
evolution of thi s threshold is not well-determined (but see 
ICoil et aLll2006al who find a similar "group scale" halo mass 
at z = 1), nor is the rate or behavior of merging within such 
halos. An a priori theoretical model for the prevalence of 
quasars in halos of this mass is therefore outside the scope 
of this paper, but remains an important topic for future work. 

Since it is also established, as discussed in § [2] that the 
characteristic mass of active BHs increases with redshift, this 
implies substantial evolution in the ratio of BH to host halo 
mass to z ~ 2. It is unclear how much of this may owe to evo- 
lution in the ratio of BH to host stellar ma ss: observational es- 
timates imply some such evolutio n (e.g.. IShields et alj|2005t 
r oo et al . 2006), b ut upper limits fro m evo- 



timates imply some 
Pen g et al.ll2006t Iw 



lution in stellar mass densities dHopkins et a l. 2006d) allow 
only a factor ~ 2 evolution by z = 2. Therefore, there might 
also be at least some increase with redshift in the characteris- 
tic ratio of stellar to halo mass. Future constraints from halo 
occupation models or galaxy clustering at high redshifts will 
be valuable in breaking this degeneracy, and potentially pro- 
vide important clues to galaxy assembly histories. 

We contrast these predictions with two extremely sim- 
ple models. In the first, quasar activity is an unbiased 
tracer of dark matte r, i.e. b(z) = 1. This does, after all, 
appear true at z = dWake et al.l 12004 iGrazian et al.1 12004 
IConstantin & Vogelevl 120061) . This is immediately strongly 
ruled out: there is an unambiguous trend that higher-redshift 
quasars are more strongly biased (as noted in essentially all 
observed quasar correlation functions). 

Next, we consider the possibility that quasars live in the 
same halos at all redshifts. This is equivalent to some "classi- 
cal" i nterpretations o f pure luminosity evolution in the QLF 
(e.g., iMathezI [1976); i.e. that quasars are cosmologically- 
long lived (alth ough other ob servations demand a lifetime 
< 10 7 yrs; e.g. Martini 2004, and references therein), and 
dim from z = 2 to the present. It is also equivalent to say- 
ing that quasars are triggered, even for a short time, in the 
same objects over time (e.g., stochastically or by some cyclic 
mechanism). In this case, the quasar lifetime can still be 
short (with a low "duty cycle" 6 ~ 10~ 3 ), although Edding- 
ton ratios must still tend to increase with redshift. Then, 
the halo bias ev olves as Equation CD , from a z = value 
b(0) ~ 1.0-1.2 dNorberg et aLlfcOoC Although this model 
is qualitatively consistent with some quasar observations, it 
is not nearly sufficient to explain the evolution of clustering 
amplitude with redshift and is ruled out at very high (> 10<r) 
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FIG. 6. — Clustering of quasars as a function of redshift (colored points in left panels, black in other panels; as in Figure[TJ, compared to different models of 
possible host populations. Upper and lower panels plot bias and comoving correlation length, respectively. Solid line inverts the comparison in FigurefT] i.e. uses 
the estimated local clustering of red galaxies to predict the quasar clustering assuming quasars are the progenitors of present ellipticals (long dashed lines show 
~ 1 <T range from uncertainties in local bias and observed bright quasar Eddington ratios). Dotted lines correspond to halos of constant mass (as labeled). Center 
panel compares this with the observed clustering of early-type galaxies (at the characteristic red galaxy M* or L* at z < 1 ; at higher redshift fo(M* ) is not longer 
well determined, so various passive galaxy surveys are considered). Right panels compare with late-type galaxies. 



sign ificance. As noted in previous studies of quasar cluster- 
ing (ICroom et al.l2005h . quasars at different redshifts must re- 
side in different parent halo populations; quasars cannot, as a 
rule, be long lived or recurring/episodic/cyclic (although this 
does not apply to very low-accretion rate activity, perhaps as- 
sociated with "radio modes"; see e.g. Hopkins, Hernquist, & 
Narayan 2005). 

Rather than a uniform population of halos at all redshifts, 
what if quasars uniformly sample observed galaxy popula- 
tions? It is, for example, easy to modify the above sce- 
nario slightly: quasars are cosmologically long-lived or uni- 
formly cyclic/episodic, but only represent the present/extant 
population of BHs (equivalently, the present population of 
spheroids). In this case, quasar correlation functions should 
uniformly trace early-type correlations at all redshifts. 

Figure [6] compares observed early-type/red galaxy cluster- 
ing as a function of redshift with that measured for quasar 
populations. At low redshifts z < 1, both mass functions 
and clustering as a function of mass/luminosity are reason- 
ably well-determined, so we plot clustering at the charac- 
teristic early-type (Schechter function) M* or L„. At higher 
redshift, caution should be used, since this characteristic 
mass/luminosity is not well determined, and so we can only 
plot clustering of observed massive red galaxies (which, 
given the observed dependence of clustering amplitude on 
mass/luminosity and color, may bias these estimates to high 
b(z) if surveys are not sufficiently deep to resolve M* or L*). 
There is also the additional possibility that the poorly known 
redshift distribution of these objects may introduce artificial 
scatter in their clustering estimates. Bearing these caveats in 
mind, the clustering of quasars and red galaxies are inconsis- 



tent at high (> 6a) significance. Quasars do not uniformly 
trace the populations of spheroids/BHs which are "in place" 
at a given redshift. Note, however, that in this comparison the 
systematic errors almost certainly dominate the formal statis- 
tical uncertainties, so the real significance may be consider- 
ably lower. 

An alternative possibility is that BH growth might uni- 
formly trace star formation. In this case, quasar cluster- 
ing should trace the star-forming galaxy population. Fig- 
ure [6] compares observed late-type/blue/star-forming galaxy 
clustering as a function of redshift with that observed for 
quasar populations. Again, at z < 1 we plot clustering at 
the characteristic M„ or L». At higher redshift we can only 
plot "combined" clustering of observed star-forming popu- 
lations (generally selected as Lyman-break galaxies); again 
caution is warranted given the known depende nce of clus- 
tering on galaxy mass/luminosity (for LBGs, see lAllen et ail 
2005). In any case, the clustering is again inconsistent at 
high (> lOcr) significance. Quasars do not uniformly trace 
star-forming galaxies . This appears to be contrar y to some 
previous claims (e.g., Adelberger & Steidel 2005a); however, 
in most cases where quasars have been seen to cluster simi- 
larly to blue galaxies, either faint AGN populations (not ~ L* 
quasars) or bright (^ L*) blue galaxies were considered. In- 
deed, quasars do cluster in a manner similar to t he bright- 
est blue galaxies observed at several redshifts (e.g.. lCoil et al.l 
l2006d: lAllen etal.ll2005l at z ~ 1 and z > 2, respectively). 
This should not be surprising; since quasars require some cold 
gas supply for their fueling, they cannot be significantly more 
clustered than the most highly clustered (most luminous) pop- 
ulation of galaxies with that cold gas. Again, this highlights 
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FIG. 7. — As Figure[6](upper left), but showing only single quasar clustering 
measurements to highlight the significance of these comparisons from any 
individual survey. The two 2dF results are not independent, but use different 
methods to derive the quasar bias. 

the fact that the real systematic issues in this comparison prob- 
ably make the significance considerably less than the formal 
~ 10 cr seen here. 

We would also like to compare quasar clustering di- 
rectly to the clustering of gas-rich (luminous) mergers. Fig- 
ure [8] attempts to do so, using available clustering measure- 
ments for likely m ajor-merger populations. At low redshifts, 
iBlake et ail d2004t) have measured the clustering of a large, 
uniformly selected sample of post-starburst (E+A or K+A) el- 
liptical galaxies in the SDSS, which from their colors, struc- 
tural properties, and fading morphological disturbances (e.g., 
lGotol l2005, and referenc es therein) are belie ved to be recent 
major merger remnants. Inf ante et alj (|2002) have also mea- 
sured the large-scale clustering of close galaxy pairs selected 
from the SPS S at low redshift. At high redshift, no such sam- 
ples exist, but Blain et aTl d2004l) have estimated the clustering 
of a moderately large sample of spectroscopically identified 
sub-millimeter galaxies at z ~ 2-3, for which the similarity 
to local ULIRGs in high star formation rates, dust content, 
line profiles, and disturbed morphologies sugges ts they are 
systems undergoing major, gas-rich mergers (e.g. JPope et alJ 
2005, 2006; Chakrabarti et al. 2006, and references therein). 
The clustering of these populations is consistent at each red- 
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FIG. 8. — As Figure|6] but comparing observed quasar clustering (colored 
points) as a function of redshift to that various populations usually associ- 
ated with galaxy mergers (black points): post-starburst (E+A/K+A) galaxies, 
close galaxy pairs, and sub-millimeter galaxies. 

shift w ith observed quasar clustering (see also Hopki ns et al. 
2007a). This is contrary to the conclusions of Blai n et al. 
(2004), but they compared their SM G clustering measu rement 
with earlier quasar clustering data dCroom et al.l l2002) below 
their median redshift z ~ 2.5. Figure [8] demonstrates that the 
dependence of quasar clustering on redshift is such that at the 
same redshifts, the two agree very well. However, given the 
very limited nature of the data, and the lack of a uniform se- 
lection criteria for ongoing or recent mergers at different red- 
shifts, we cannot draw any strong conclusions from the direct 
merger-quasar clustering comparison alone. 

A lthough quasars do not ap pear to trace star-forming galax- 
ies, lAdelberger et alj (I20051 and references therein) have 
shown that the star formation rates, clustering properties, 
and number densities of high-redshift LBGs suggest they are 
the progenitors of present-day ellipticals. To the extent that 
quasars are also the progenitors of ellipticals (but with a larger 
clustering amplitude at a given redshift compared to LBGs), 
this suggests a crude "straw-man" outline of an evolutionary 
sequence with time, from LBG to quasar to remnant elliptical 
galaxy. Knowing how the clustering properties of halos host- 
ing LBGs with a given observed bias at some redshift zlbg 
will subsequently evolve, we can determine the redshift zq 
at which this matches observed quasar clustering properties. 
This offset, if LBGs and quasars are indeed subsequent pro- 
genitor "phases" in the sequence of evolution to present day 
ellipticals, defines the "duration" of the LBG "phase" or time 
between LBG and quasar "stages." 

Figure [9] considers this in detail. We show the observed 
clustering of quasars and LBGs from Figure|6l with curves il- 
lustrating the subsequent clustering evolution of the LBG host 
halos observed at z = 1, z = 2, and z = 3. These correspond to 
the characteristic observed quasar clustering at z = 0.4, z = 1 .0, 
and z = 1.3, respectively. Thus, halos of the characteristic 
LBG host halo mass at z = 3 will grow to the characteristic 
quasar host mass at z = 1.3, and so on. We also show the 
physical time corresponding to this offset, calculated from the 
observed LBG clustering at various redshifts and the best-fit 
estimate of the LBG host mass ~4x 10" /r'M Q , and this 
time divided by the Hubble time (age of the Universe) at the 
"quasar epoch" zq- Interestingly, this implies that objects 
characteristically spend ~ 3-4Gyr (~ l/2t H at the redshifts 
of interest) in the "LBG phase." This may reflect the time for 
dark matter halos to grow from the characteristic LBG mass, 
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at which star formation and the conversion of mass t o light 
appears to be most efficient (e.g., Whit e & Reesl [l978) to the 
typical quasar host mass; but it is also possible that associated 
physical processes related to quasar fueling or the termina- 
tion of star formation set this timescale. If quasars are trig- 
gered in major mergers, this rather large time offset (as op- 
posed to the typical ~ lOOMyr delay b etween starburst and 
quasar in major merger simulations, Di M atteo et al.l 120051) 
implies that LBGs are themselves not primarily driven in ma- 
jor mergers. A similar conclusion was recently reached by 
lLaw et all d2006l) from direct analysis of LBG morphologies 
at z ~ 2-3. This conclusio n and the LBG clustering in Fig- 
ure [9] ( Wechsle r et al.l 1200 ll) are broadl y consistent with the 
expectations of semi-analytic models (Some rville et al.ll2001l) 
which argue LBGs are driven largely by "collisional" minor 
merging. 

We can also use this approach to determine the time be- 
tween the "quasar" and red/elliptical phases in this evolution- 
ary sequence. Figure [TOl shows this, in the style of Figure |9l 
where the redshift shown in the middle and lower panels refers 
to the redshift of the observed quasar population, and the time 
to the delay at which their evolved clustering matches that 
measured for the red galaxy population. Note that the contin- 
uous curves calculated in the middle and lower panels assume 
the red galaxy clustering is well-fitted by the plotted (upper 
panel) constant halo mass ~ 1.6 x 10 13 h~ l M@ curve; this is, 
in fact, not a very good approximation at low redshifts, hence 
these curves diverge below z ~ 1 -2 from the times calculated 
from the actual red galaxy clustering measurements. 

In the lower panels, we also plot the time for "burst- 
quenched" star_formation history models adapted from 
lHarker et al.l (|2006) to redden to a typical constant "red 
galaxy" threshold rest-frame color U—B > 0.35. These model 
star formation histories assume a constant star formation rate 
until 1 Gyr before the "quasar epoch," then a factor 5 enhanced 
star formation rate for this 1 Gyr, at which point star forma- 
tion ceases. Essentially, this yields a useful toy model for 
"quenching," if indeed the triggering of quasars is associated 
with the formation of ellipticals or termination of star for- 
mation (the pre-quenching enhancement being an approxima- 
tion to , e.g., merge r -indu ced star formation enhancements), 
which [Harker et al.l (120061) demonstrate yields a reasonable 
approximation to the observed mean color, number density, 
and Balmer H<5p absorption strength evolution of red galax- 
ies. The predicted time for such quenched star formation his- 
tories to redden to typical red galaxy colors agrees well with 
the time estimated from clustering here at all redshifts; i.e. 
the color and halo mass evolution of these systems are consis- 
tent with reasonable star formation histories in which quasar 
activity is associated with "quenching" or the termination of 
star formation. 

We have estimated the time offsets in Figures [9] & [10] 
from a direct comparison of the observed clustering. In- 
stead, one might imagine adopting the implied halo mass 
(~4x 10 h Mq) at the "star-forming" phase and using ex- 
tended Press-Schechter (EPS) theory to calculate the average 
time for a typical progenitor halo of this mass at each ob- 
served redshift to grow to the implied quasar host halo mass 
(~4x lO 12 /; -1 Mq). We discuss this in greater detail in §|5J 
and show that it has no effect on our conclusions. For the pur- 
poses here, adopting this methodology (specifically, calculat- 
ing the evolution of the "main branch " proge nitor halo mass 
with redshift following iNeistein et al.l (120061) in our adopted 
cosmology) systematically increases the inferred time delays 
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FIG. 9. — Upper: Clustering of quasars and star-forming galaxies, as in 
Figure|6] Solid green lines show the subsequent evolution of the clustering of 
the star-forming galaxy halos from z = 1, 2, 3. Middle: Time delay between 
the star-forming or LBG phase and "quasar" epoch, defined as the time after 
the observed redshift of each LBG population at which its evolved clustering 
will match that of the obs erved quasar population. Points as in upper panel, 
with data from Ouchi et al. (2004, cyan inverted triangles) added at z > 4. 
Dashed and dotted lines show time for halos of mass ~4x 10 1 1 hr l Mq (the 
typical LBG host mass) at each redshift to reach ~4x 10 12 h" 1 Mq (dotted) 
or the (weakly redshift-dependent) halo mass defined by our best-fit trend in 
the upper panel (solid). Lower: As the middle panel, but the time shown is as 
a fraction of the Hubble time at the "quasar" epoch. 



(points) in Figure |9]by ~ 1 - 2 Gyr and those in Figure [TOlbv 
~ 0.5 - 1 Gyr, but does not significantly change the plotted 
trends or comparisons. 

So, this leaves us with the following suggested empirical 
picture of galaxy evolution. Galaxies form and experience a 
typical "star forming" or LBG epoch, with maximal efficiency 
around a characteristic halo mass of a few - IO'^-'Mq. 
Growth continues, presumably via normal accretion, minor 
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FIG. 10. — As Figure[9] but instead showing the time from the "quasar 
phase" to the "red galaxy" phase implied by the observed clustering of both 
populations. Green long-dashed lines in the middle and lower panels show 
the t ime required for th e "burst-quenched" star formation history models 
from Harker et al. 1 2006) (which yield a good empirical approximation to the 
buildup and mean color evolution of red galaxies) to redden to a threshold 
U-B> 0.35. 

mergers, and star formation, for roughly half a Hubble time, 
until systems have growth to a characteristic halo mass ~ 
4 x 10 l2 h~ l Mq. At this point, some mechanism (for exam- 
ple, a major merger, as this may be the characteristic scale 
at which the host halo grows large enough to host multiple 
"large" star-forming systems) triggers a short-lived "quasar" 
phase, drives a morphological transformation from disk to 
spheroid, and terminates star formation. About ~ 1 — 2 Gyr 
after this, the host halos have grown to ~ 1O i3 /i -1 M and the 
spheroids have reddened sufficiently to join the typical "red 
galaxy" population. They then passively evolve (although 
they may experience some gas-poor or "dry" mergers) to z = 0, 
satisfying observed correlations between BH and spheroid 



properties. Although individual BHs can, in principle, gain 
signifi cant mass from "dry" mergers (see, e.g. iMalbon et alJ 
2006), this cannot (by definition) add to the total mass budget 
of BHs, which must be built up via accretion. Note that this is 
only a rough conception outline of an "average" across popu- 
lations and should not be taken too literally. Different systems 
will undergo these processes at different times, and (possibly) 
via different mechanisms. Still, this provides a potentially 
useful framework in which to interpret these observations. 

5. AGE-MASS RELATIONS AND CLUSTERING 

In Figure Q~T] we compare the mean age of BHs of a 
given z = mass with that of the stellar population of their 
hosts. At a given redshift, the characteristic QLF luminos- 
ity L*(z) and corresponding characteristic "active" mass Mbh 
from Figure [2] define the epoch of growth of BHs of that 
mass. The typical "age" of BHs of that mass will be the 
time since this epoch. In detail, Equation © relates the 
observed QLF to the rate at which BHs of a given relic 
mass are formed as a functi on of redshift. We ad opt the 
fits to this equation given in Hopkins et al. (200 7b(). which 
use the model quasar lightcurves determined in Hopkins et al. 
(2006a,b) to calculate the time-averaged rate of formation of 
individual BHs as a function of mass and redshift, to esti- 
mate the median age (peak in rate of formation/creation of 
such BHs as a function of time) and 25-75% interquartile 
range in "formation times." This introduces some model de- 
pendence, but as discussed in § [2] a similar result is obtained 
using very different me thodologies, including purely empiri- 
cal, simplified models (lYu & T remaine 2002c lMerlonill206"4l 
iMarconi et~ai1l2004 [Shankar et alJl2004l) . or di rect calcula- 
tion from observed Eddington ratio distributions (jVeste rg aard 
12004 iMcLure & Dunlopll2004 iKollmeier et al.ll2005k 

In any case, we recover the well-known trend that the more 
massive BHs are formed at characteristically earlier times 
dSalucci et alj|1999t lYu & Tremaind 120021: lUeda et alj | 2003b 
Heckman et all 12004 lHasinger et al] 12005b iMerlonil 12004 
Marconi etall 120 04: Sha nkar et alj 12004 iMcLure & Dunlopl 
2004 lKollmeie r~t alll2005l) . This is not surprising, as most 
massive BHs must be in place by z ~ 2 to power the brightest 
quasars, and these objects are generally "dead" by low red- 
shift (with lower-mas s objects dominating the local QLF, e.g. 
iHeckman et aT]|2004l) . 

Given a BH mass, we can compare with the observation- 
ally determined age of the typical host galaxy (with Mbh = 
/LtMgai). First, we consider early-type hosts, specifically the 
stellar ages of host spheroids of BHs at each mass Mbh- The 
mean ages (and dispersion about that mean) of ellipticals as 
a function of stellar mass have been estimated in a num- 
ber o f studies, recently for example by Gallazzi et al] d2005l 
2006), who fit SDSS spectra (line indices) and photometry 
for ~ 175,000 local galaxies to various realistic star forma- 
tion histories, including a mix of continuous and/or starburst- 
ing histories while allowing mass, total metallicity, and abun- 
dance ratios to freely vary. They quote r-band light-weighted 
ages, which for our purposes are effectively equivalent to the 
ages determined by fitting a single stellar population (SSP) or 
"single burst" model to observed spectra, and indeed agree 
very well with best-fit SSP ages from simila r SDSS sam- 
ples (IClemens et al]|2006l: iBernardi et al . 2006) a nd previous 
studies ( J0rgensen 1999; Trager et al. 2000; Kuntschner et al. 
1 2001 1 : [Caldwell et al 2003; B ernardi et al. 2005; Nela n et all 
120051: iThomas et al.ll2005b Icollobert et al.ll2006l for a review 
see Renzini 2006) as a function of elliptical stellar mass. A 
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FIG. 11. — Left: Upper panel shows the mean z = age (lookback time to the mean "formation epoch") of BHs as a function of mass (black solid line, dashed 
lines show 25% — 75% quartile r anges), compared to the stellar populati on age of their hosts. Ages of spheroids as a function of ma ss (with Mbh = M^gal) 
are shown (colored points) from Nelan et al. (2005, NFPS; red squares), Gallazzi et al. (2006, SDSS; blue stars), Thomas et al. (2005, orange circles; "field" 
subsample). Errors show the dispersion in ages at a given mass. Lower panel uses this age to predict quasar clustering as a function of redshift; i.e. assuming the 
"quasar epoch" of spheroids of a given mass is associated with the termination of star f ormation (black lines, as labele d; colored points sho w observed quasar 
clustering as in Figure[6}. Center: Same, for ages of host disks; ages from r-model fits of Bell & de Jong 1 2000, red) and Gavazzi et al. (2002, orange) (the offset 
between them owes to the choice of initial time in the r-model). Solid lines asumme Mbh <x Mb u i ge , dotted Mbh <x iM&A. +Mb u i ge ). Dashed lines re-calculate 
the age for a single-burst SFH. Lower panel is as lower left, assuming quasar acti vity is associated with the star formation epoch (as labeled). Right: Blue (solid) 
line shows the "all progenitor" age (DM "downsizing" from Neistein et al. 2006), red (dashed) the age of the main progenitor halo, and green (dotted) the time 
when halo crossed the "quenching" mass from Dekel & Birnboim 1 2006). Lower panel as lower left, assuming quasar age is equal to the halo age as labeled. 

range in host masses at a given time, the agreement is rea- 
sonable. Including the dispersion in ages, i.e. modeling the 
age distribution at each Mbh as a Gaussian with the observed 
scatter, improves the agreement and yields a nearly identical 
prediction of bias as a function of redshift to that in Figure [6] 
(solid black line). 

We can of course repeat these exercises for other pos- 
sible "host" populations. We next consider correlations 
with the star formation histories of late-type or star-forming 
host galaxies - i.e. the possibility that quasar activity is 
generically associated with star formation. The observed 
star formation histories are similarly estimated, generally 
by fitting to exponentially declining models (r-models; star 
formation rate M oc exp[-(f-f,)/r] since an initial cos- 
mic time ti). Specifically, we consider t he fits of late- 
type a ges as a function of stell ar mass from Bell & de Jons 
(2000) and Gavazzi et al. ( 200 2J) (consisten t with Janse n et alJ 
| 2000t iBell et alj l2000t iBoselli et al] 1200 U Ik auffman n et all 
2003at iPerez-Gonzalez et alJ 120031: fBrinchmann et alJl2004t 
iMacArthur et al.1 12004 iGallazzi et all l2005h . The mass- 
weighted age is c alculated fro m the model SFR (see 
[ Bell & de JongT2 000. Equation 3). IBell & de Jongl d2000t) and 
iGavazzi et a l. (2002) technically quote the age and metallic - 
ity as a function of K and H band absolute magnitudes, re- 



similar res ult is also obtained independen tly by iTreu et ail 
d2005l) and Idi Serego Alighieri et all d2006l) from studies of 
the the fundamental plane evolution of early-type galaxies. 
Note that the error bars shown are the measured dispersions 
in the population about the mean age at a given mass, not the 
uncertainties in the mean ages themselves (which are smaller; 
~ 0.2 Gyr statistical, ~ 1 Gyr systematic; see iNelan et alj 
120051) . The agreement between BH and host stellar ages is 
good at all masses; both the trend and dispersion (interquar- 
tile or ±lcr range) about it are similar (\ 2 jv ~ 8/17 for a 
direct comparison). 

If the age of its stellar populations is indicative of when the 
"quasar epoch" occurred in a given host, then, without mak- 
ing any assumptions about the masses of black holes or quasar 
Eddington ratios, we can use the mean age of stellar popula- 
tions to predict quasar clustering. In this scenario, ellipticals 
of mass M, with mean age fh OS t, would represent the popula- 
tion "lighting up" as quasars at a lookback time of fhost, and so 
the quasar bias at that lookback time should be the local bias 
of ellipticals of mass M (Equation|7]l, evolved to the appropri- 
ate lookback time fh os t with Equation ([]]). FigurefTTIcompares 
this expectation with the observed quasar bias as a function of 
redshift. Despite the very simple nature of this model, which 
ignores both the range of ages at a given M and, similarly, the 
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spectively, but given their quoted best-fit stellar population 
models at each luminosity, it is straightforward to calculate 
the corresponding mass-to-light ratios {M/Lk and M/Lh) and 
convert the observed luminosities to total stellar masses. To 
convert to a corresponding BH mass, we consider first a uni- 
form application of the local BH-host mass relation, i.e. as- 
suming BH mass is correlated with total stellar mass, and sec- 
ond determining the mean bulge-to-disk ratio for a given total 
late-type stellar mass or l uminosity ( seelFukugita et alj[l998l; 
I After & Richstond 120021: iHunt et al] 12004 for the appropri- 
ate mean B/T for different masses/luminosities) and assum- 
ing BH mass is correlated with the bulge mass only. Because 
the trend of age as a function of stellar mass is weak, consid- 
ering the total mass or bulge mass makes little difference in 
our comparison, and we subsequently consider the observa- 
tionally preferred correlation between BH and bulge mass. 

The mean age of a given population derived from different 
model star formation histories will, of course, be weighted 
differently. To show the systematic effects of such a choice, 
we roughly estimate the equivalent age from a single burst or 
SSP model. We calculate the z = observed (B-V) color at 
each mass from the mean best-fit r models, and then calcu- 
late the corresponding age for the same {B-V) of a single 
burst model (of th e same metallicity as a fu nction of mass) 
from t he models of lBruzual & Chariot! d2003) with a Salpeter 
d 19551) IMF. Although most of the observations above find 
this SSP approximation is not good for star-forming galaxies, 
it illustrates an important point. The SSP ages are weighted 
towards the youngest, bluest stellar populations - essentially 
functioning as an indicator of the most recent epoch of signif- 
icant star formation, and are therefore quite young (< 2Gyr; 
similar to the typical time since recent low-level starbursting 
activity found in late-types with the mor e realistic star forma- 
tion histories in Kauffmann et al. 2003a). However, the trend 
as a function of mass is unchanged and the overall agreement 
is worse. Therefore, while the systematic effects here are sub- 
stantially larger than the measurement errors in mean age as 
a function of mass, they cannot remedy the poor agreement 
with the ages of BH populations. 

We again use this age as a function of total/bulge mass, and 
the observed clustering of late-type galaxies from Figure Q] at 
z = 0, to estimate what the quasar bias as a function of red- 
shift should be, if these systems were the hosts of quasars 
and their quasar epoch were associated with the age of their 
stellar populations. The predictions are inconsistent with the 
observations at high significance (> 4.5c), regardless of the 
exact age adopted (r-model or SSP). In fact, the predicted 
clustering as a function of redshift is highly unphysical (ow- 
ing to the fact that there is relatively little difference in ages, 
but strong difference in clustering amplitudes from the least 
to most massive disks). Ultimately, this demonstrates that 
the hypothesis that quasar activity generically traces star for- 
mation is unphysical. This is also supported by the fact that 
the integrated global star formation rate and quasar luminos- 
ity density ev olve in a similar, but not identical manner from 
z ~ 0-6 (e.g jMerloni et al.ll2004l) . 

We next consider the the possibility that quasar activity 
traces pure dark matter assembly processes - i.e. that the 
buildup of BHs in quasar phases purely traces the formation 
of their h ost halos. Given the loca l BH-host stellar mass rela- 
tion from Mar coni & Hunj d2003l) . and the typical halo mass 
as a function of early-ty pe hosted galaxy mass cali brated from 
weak lensing studies by Mandelbau m et al] (120061) . we obtain 
the mean host halo mass as a function of BH mass (mean 



A^haio ~ 4 x 10 4 Mbh; although the relation is weakly non- 
linear). For our adopted cosmology, we then calculate the 
mean age (defined as the time at which half the present mass 
is assembled) of the main progenitor halo for z = halos of 
this mass. Error bars are take n from an ensemble o f random 
EPS merger trees following Nei stein et alj (120061) . Know- 
ing the mass of a h alo at a given redsh ift, we calculate its 
clustering following [Mo & Wh ite (1996) as in § [4] and use 
this combined with the mean ages to estimate the expected 
quasar clustering as a function of redshift if quasars were as- 
sociated with this formation/assembly of the main progenitor 
halo. Although the exact age will depend on cosmology and 
the adopted "threshold" mass fraction at which we define halo 
"age," the result is the same, namely we recover the well- 
known hierarchical trend in which the most massive objects 
are younges t, in contradiction wit h quasar/BH ages. 

However, Nei stein et al] (120061) have pointed out that the 
mean assembly time, considering all progenitor halos, can 
exhibit so-called "downsizing" behavior. We therefore fol- 
low their calculation of the mean age of all progenitors as 
a function of z = halo (and corresponding BH) mass, and 
also use this to estimate quasar clustering as a function of 
redshift. Again, although the systematic normalization de- 
pends somewhat on our definitions, it is clear that the recov- 
ered "downsizing" trend is, as the authors note, a subtle ef- 
fect, and not nearly strong enough (inconsistent at > 10<r) 
to explain the downsizing of BH growth. Again, the abso- 
lute value of the age obtained can be systematically shifted by 
changing our definition of halo "assembly time," but the trend 
is not changed and significance of the disagreement with BH 
formation times is still high. 

Certain feedback-regulated models predict that black hole 
mass should be c orrelated with halo circular velocity (a s 
M BH oc v 5 c or oc v 4 c : iSilk & Reel[l998t IWvithe & Loebll2003l) . 
rather than halo mass. To consider this, we have re-calculated 
the "all progenitor" and "main progenitor" ages, but instead 
adopted the time at which the appropriate power of the circu- 
lar velocity (v^ or v 4 ) reaches half the z = value. Because, 
for a given halo mass, v c is larger at higher redshift, this sys- 
tematically shifts both ages to higher values, but the trends 
are similar. In each case the resulting ages disagree with the 
quasar/BH ages at even higher sig nificance. 

Alternatively, some models (jBirnboim & Dekel 120031; 
iBinnevI 120041 iKeres et alj 120051 iDekel & Birnboiml [2006b 
suggest that a qualitative change in halo properties occurs at 
a characteristic mass, above which gas is shock-heated and 
cannot cool efficiently, forming a quasi-static "hot accretion" 
mode in which quasar feedback can act efficiently. Although 
there is no necessary reason why quasar activity should be 
triggered by such a transition, its posited association with 
quasar feedback leads us to consider this possibility. 

Knowing the mean z = halo mass for a given Mbh, we 
plot the age at which the main progenitor halo mass surpassed 
the crit ical "quenching" mass def ined as a function of red- 
shift in Dekel & Birnboim (2006). Since this amounts to a 
nearly constant characteristic halo mass ~ x 10 11 - 10 12 M Q , 
the expected quasar clustering as a function of redshift is not 
unreasonable (see Figure [6l there is a systematic offset, but 
this is sensitive to the adopted cosmology). However, this 
model actually predicts too steep a trend of age with mass 
(inconsistent at > 6a). The ages of the most massive sys- 
tems are reasonable (which, in com parison to the ages of el- 
lipticals, has been widely discussed; IDekel & Birnboiml2006l: 
ICroton et al.ll2006t iDe Lucia et aIB 2006). but the host halo of 
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FIG. 12. — Observed clustering of quasars, compared to that inferred from 
their z = early-type hosts if the termination of st ar fo rmation is coincident 
with quasar activity (as lower left panel in Figure [TTJ. Left: Our standard 
methodology is used to empirically evolve the clustering of local systems 
(black points) to the redshifts shown. Right: Instead, using the full EPS 
formalism and estimated b(M^\ , z) to evolve the clustering of local systems. 
Differences owing to the choice of methodology are small at the halo masses 
of interest. 

a typical ~ 10 8 M Q BH (i.e. SchechterM*) crosses the thresh- 
old halo mass a mere ~ 5 Gyr ago, predicting, in this simple 
model, that these BHs should have been the characteristic ac- 
tive systems at z ~ 0.5 instead of the observed z ~ 1.2—1.5. 
At lower masses, the mean z = halos are only just at, or are 
still below, this critical halo mass. Given the scatter in the BH- 
host mass relations, there will be some BHs of these smaller 
masses living in larger halos which have already crossed the 
"quenching" mass, but the age distribution will still be one- 
sided and weighted to very young ages. To match the ob- 
served age trend, there must, in short, be some process which 
can trigger quasar activity at other halo masses before they 
cross the "quenching" threshold. 

Finally, we note that in evolving the clustering of local sys- 
tems "up" in redshift in the lower panels of Figure QT| there 
might be some ambiguity (if, for example, a given z = halo 
is assembled from many progenitor high-redshift halos with 
significantly different properties). The simple evolution pre- 
dicted by Equation (Q~|) is derived from pure gravitational mo- 
tions, and therefore as applied moving "backwards" in time 
represents an ef fective "m ean" bias of the progenitors of the 
z = system (see Fry 1996). To the extent, however, that there 
is a dominant progenitor halo at a given redshift and many 
smaller halos which will be accreted by the "main" halo, it 
is the properties of the main progenitor which are of interest 
here. 

We therefore consider a completely independent approach 
to empirically compare the clustering measurements shown, 
which attempts to capture these subtleties. Given a z = popu- 
lation, we can estimate its characte ristic host halo mass eithe r 
directly from the measurements of lMandelbaum et aH (2006), 
or indirectly by matching the observed bias (with bias as a 
function of halo mass calculat ed fo r the adopted cosm ology 
following iMo& Whitd d 19961) and ISheth et al.1 ( 12001 ) as in 
§ HI. Following iNeistein et alJ (120061) . we then calculate the 
mass of the main progenitor halos of this z = mass, as a 
function of redshift (i.e. the highest-mass "branch" of the EPS 
merger tree at each redshift). At the redshift of interest (e.g. 
appropriate lookback time, for the comparisons in the lower 
panels of Figure [TTJ, we then calculate the expected bias for 
halos of this main progenitor mass. 

Figure [12] reproduces the lower-left comparison in Fig- 
ure [TTJ (the expected clustering of elliptical progenitors at 
the times determined by their stellar population ages), us- 



ing both our previously adopted methodology and this re- 
vised estimation. The latter method has the advantage, as 
noted above, of accounting for the difference between the 
main progenitor and smaller, accreted systems. The ap- 
proach, however, suffers from certain inherent ambiguities 
in Press-Schechter theory. For example, the calculated evo- 
lution is not necessarily time-reversible, and the clustering 
properties are assumed to be a function of halo mass alone, 
which recent high-resolu tion numerical simulations suggest 
may not be correct (e.gjGao & Whitell2006l : |G. Harker et alJ 
l200aiWechsleretal. 2006). In particular, if quasars are trig- 
gered in mergers (i.e. have particularly recent halo assem- 
bly times for their post-merger halo masses), then they may 
represent especially biased regions of the density distribu- 
tion. Unfortunately, it is not clear how to treat this in de- 
tail, as there remains considerable disagreement in the liter- 
ature as t o whether or not a significant "merger bias" exists 
(see, e.g. [K auffman n & Haehneltll2002t iPercival et ail 120031: 
iFurlanetto & Kamionkowski 2006). Furthermore the dis- 
tinction between galaxy-galaxy and halo-halo mergers (with 
the considerably longer timescale for most galaxy mergers) 
means that it is not even clear whether or not, after the galaxy 
merger, there would be a significant age bias. In any case, 
most studies su ggest the effect is quite sma ll: using the fitting 
formulae from Wechsler etin (|200l[2006), we find that even 
in extreme cases (e.g. a M 3> M v ; r halo merging at z = as 
opposed to an "average" assembly redshift z/ « 6) the result 
is that the "standard" EPS formalism underestimates the bias 
by ks 30%. For the estimated quasar host halo masses and 
redshifts of interest here, the maximal effect is < 10% at all 
Z = — 3, much smaller than other s ystematic effe c ts we have 
considered. This is consistent with iGao & W hite (2006)) and 
ICroton et al.1 (120071) who find that "assembly bias" is only im- 
portant (beyond the 10% level) for the most extreme halos or 
galaxies in their simulations. 

In practice, Figure[T2]demonstrates that, for the halo masses 
of interest here, the two methods yield very similar results. 
This is reassuring, and owes to the fact that the differences 
from the choice of methodology discussed above are impor- 
tant only at very high or very low halo masses, where for 
example the clustering of small halos which are destined to 
be accreted as substructure in clusters (> 10 l5 /i _1 M Q ) will 
be very different from the clustering of similar-mass halos in 
field or void environments. Alternatively, one can think of 
the EPS approach as attempting to account for the possibil- 
ity that bias is a non-monotonic functi on of mass (e.g. ris - 
ing galaxy bias at very low luminosities, Norberg et al. 2002), 
which FigureUJdemonstrates is important only at masses well 
below those of interest here. To the extent that any "merger 
bias" is permanent or long-lived (as expected if the excess 
clustering is correlated with halo concentration or formation 
time), our "standard" methodology should account for it, as 
we simply evolve the clustering of the present hosts of quasar 
"relics" to earlier times according to gravitational motions. 
That the different seen is small provides a further reassur- 
ance that the effects of "merger bias" are probably not dra- 
matic. Ultimately, we have re-calculated all the results herein 
adopting the more sophisticated (but more model-dependent 
and potentially more uncertain) EPS approach, and find that 
it marginally improves the significance of our conclusions but 
leaves them qualitatively unchanged. 

6. IMPLICATIONS FOR HIGH-REDSHIFT CLUSTERING 
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At z > 2-3, comparing quasar and early-type clustering be- 
comes more ambiguous. Above z ~ 2, the QLF "turns over", 
and the density of bright quasars declines. Specifically, it 
appears that the char acteristic quasar luminosity L» declines 
(Hopk ins et al.ll2007bl and references therein), at least from 
Z ~ 2 to z ~ 4.5 above which the "break" L* can no longer 
be determined. One possible interpretation of this is an ex- 
tension of our analysis for z < 2; i.e. one could assume that 
each "quasar" episode here signals the end of a BH's growth, 
which will evolve passively to z = 0. At z = 0, the tightness of 
the local BH-host relations means the hosts must have the ap- 
propriate mass and lie within the appropriate halos, to within 
a factor ~ 2 of the observed scatter. Therefore, we can adopt 
the same approach as in §|4]to use the local observed cluster- 
ing as a function of host properties to evolve back in time and 
predict quasar clustering as a function of redshift. 

Figure [13] shows the bias and correlation length predicted 
by this approach, an extension of the model (Equation we 
have considered at z < 2. Figure [14] also shows the typical 
host halo mass corresponding to the predicted clustering as 
a function of redshift (for our adopted cosmology); in this 
simple extension of the z < 3 case, the observed decline in 
the QLF L* traces a decline in the characteristic (although not 
most massive) quasar-hosting halo mass. 

However, at high redshifts, flux limits may severely bias 
clustering measurements. Although at L < L*, quasar clus- 
tering does not strongly depend on the quasar luminosity (see 
§ [3j, implying a well-d efined characteris tic active mass which 
we can adopt (see also lLidz et al.l l2006). this is not necessar- 
ily true for L > L*. In fact, Figure [3] shows (and o bservations 
may begin to see, e.g. Porciani & Norberg 2006) a steepen- 
ing of bias versus luminosity at L > L», reflecti ng the uni- 
formly high observed quasar Eddington ratios dVestereaardl 
l200ltlMcLure & Dunlopl200l lKollrneierltlTl2005h at high 
luminosities, which imply the bright end of the QLF (L ^> L*) 
becomes predominantly a sequence in active BH mass. To the 
extent that BH mass traces host mass, then, these systems re- 
side in more massive hosts and will be more strongly biased. 

In order to estimate how this will change the observed clus- 
tering, we roughly approximate this effect as follows. For a 
given flux limit at a given redshift there is a reasonably well- 
defined survey depth, to a minimum luminosity L m [ n . If this 
is sufficiently deep to resolve the QLF "break," i.e. L m [ n < L», 
then the weak observed dependence of clustering on luminos- 
ity means the observed clustering will trace that characteristic 
of ~ L» quasars - corresponding to characteristic Mbh active 
BHs and M ga i w Mbh//** hosts (our fiducial model, and the 
case for all observations plotted in Figure [T3l , However, if 
the flux limit or redshift is sufficiently high such that > 
L*, then the survey will not sample these characteristic host 
masses. In this case, we consider the bias as a function of lu- 
minosit y plotted in Figure[3] from the models of Hopki nset al.l 
d2006bl) ; lLidz et al.1 d2006). evaluated at L min at the given red- 
shift. Qualitatively, for the nearly constant Eddington ratios 
observed at L > L», L oc Mbh <x M ga i, we expect L m ; n > L* to 
correspond to an approximate minimum observed host mass, 
M mia ~ M ga i(L*) x (L rain /L») ~M Edd (L min )//i. Since the QLF 
slope is steep at L > L», objects near L m ; n or M m ; n will dom- 
inate the observed sample, and so this amounts to calculat- 
ing the clustering for this mass, instead of M ga i(L*), at the 
given redshift. We caution that this is a rough approximation 
to more realistic selection effects, but should give us some 
idea how flux limits will bias the observed clustering. 

We consider several representative flux limits, in observed- 



frame /-band, typical of optical quasar surveys (e.g. the 
SDSS), in addition to the case with effectively infinite depth 
{mi < 30). We calculate L* at the appropriate rest-frame wave- 
lengths as a function of redshift using the fits to L* from 
Hopki ns~et al.l (l2007bl) . spanning z ~ - 6 and spanning the 
relevant rest-frame wavelength intervals. At the limits of most 
current optical surveys, ;«,■ < 20.2, the QLF break L* is only 
marginally resolved at z ~ 2-3, and so above this redshift sur- 
veys are systematically biased to more massive L > L» BHs 
and higher clustering amplitudes. However, a relatively mod- 
est improvement in depth to m, < 22 would allow unbiased 
clustering estimates to be extended to z ~ 4. 

We have so far assumed BHs effectively "shut down" af- 
ter their quasar epoch - i.e. "efficient feedback" even at 
high redshifts. Although the various observations discussed 
above (Eddington ratio distributions, quasar host measure- 
ments, HOD models, black hole mass functions, and our clus- 
tering comparison) demand this be true at z < 2 - 3, there 
are no such constraints at z > 3. In other words, it is pos- 
sible that the increase in the QLF from z ~ 6 to z ~ 3 traces 
the growth of the same populations of BHs, not the subse- 
quent triggering and "shutdown" of different populations. If 
BHs at z ~ 6 continue to grow to z ~ 2 - 3 before "shutting 
down," they must live in more massive z = host galaxies 
(to preserve the tight observed BH-host mass relation), and 
thus should have stronger clustering amplitudes. We therefore 
consider two representative simple models which bracket the 
range of possibilities for this growth and present simple tests 
for future clustering measurements to break the degeneracy 
between these models. 

First, we assume that quasars grow with the QLF to z ~ 2 
before "shutting down" (i.e. "inefficient feedback"). In such 
a case, z ~ 6 quasars grow either continuously or episodi- 
cally with their host systems until the epoch where "down- 
sizing" begins, and the QLF at all redshifts z > 2 represents 
the same systems building up hierarchically. The z = relic 
masses (and therefore z = characteristic host masses, from 
which we calculate the "parent" halo clustering as a func- 
tion of redshift) are then the same at all z > 2. This is 
also equivalent to a "pur e density evolu tion" model for the 
high redshift QLF, as in iFan et al.l (1200 ll) . in which the QLF 
break luminosity remains constant above z ~ 2-3 while 
the number density/normalization uniformly declines. Such 
a model is marginal ly disfavored by current measurements 
(Hopk ins et al.l l2007b), but constraints on L* at high redshifts 
are sufficiently weak that it remains a possibility. 

The clustering as a function of redshift in this model be- 
haves very differently from the previous model at high red- 
shifts. If objects cannot grow after their quasar epoch even 
at high redshifts, then the subsequent decline of the QLF L* 
traces a decline in characteristic active masses, and the bias 
of active systems "turns over"; however, if all grow to the 
characteristic at z ~ 2, then these high redshift systems all 
represent similar z = masses by the time they "shut down," 
and must be increasingly biased at higher redshifts. 

Next, we consider a "maximal" growth model, in which 
we assume not only that the buildup of the QLF represents 
the continued growth of BHs until z ~ 2, but also that this 
growth proportionally tracks the typical growth of dark matter 
halos over this redshift range. We very crudely estimate the 
"typical growth" with the growth of an average high-redshift 
quasa r "host" halo. Based on their space density and BH 
mass, IFan et alj d2001l 120031) estimate that typical z ~ 5 - 6 
SDSS quasars represent ~ 6er overdensities. We therefore as- 
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FIG. 13. — Using the model in Figures [TJ& [6] to predict the observable clustering of quasars at high redshifts z 3. Observed clustering is shown (points) as 
in FigurefJJ Different lines show the effect of different observed-frame /-band flux limits, as labeled (note i = 20.2 corresponds to the SDSS DR3 completeness 
limit; [Richards et al. 2006a). Left panels assume efficient feedback at high redshifts; i.e. that BH growth "shuts down" after each quasar episode. Center panels 
assume all z, > 2 BHs grow with the observed QLF to the characteristic peak luminosities at z, ~ 2, then shut down ("inefficient feedback"). Right panels assume 
quasar growth tracks host halo growth, even after a quasar episode, until z. = 2 ("maximal growth"). Future observations at z ~ 4 with moderately improved flux 
limits m, < 22 should be able to break the degeneracies between these models. 

which growth "shuts down." This then yields the z = BH 
mass, corresponding host mass, and evolved clustering. Note 
that, although similar, this is not the same as assuming quasars 
track 6a overdensities at z > 2, because to the extent that the 
QLF L* does not grow by the same proportionality, this model 
effectively allows "new" or different BHs/host halos to dom- 
inate the QLF at different redshifts. It simply mandates that 
they all grow at this rapid rate. For example, an observed z ~ 6 
BH of ~ 1O 8 M is assumed to reach a mass of 2 x 1O 9 M at 
z = 2 (and then shuts down, so that this is also the mass at 
z = 0), and a ~ 10 8 M Q observed quasar at z = 4 will grow to 
~5x 10 8 M Q . The choice of rate is arbitrary, we choose it 
as a reasonable upper limit. In any case, the predicted evo- 
lution of the bias as a function of redshift is extremely steep, 
so the exact values will be very sensitive to the growth model 
and adopted cosmology. The point we wish to illustrate is that 
this model generically predicts a steep bias evolution at z > 2, 
which regardless of the details will be distinguishable if future 
quasar clustering measurements at z = 3 - 4 can be extended 
to a depth of m, < 22. 

Note that extending the depth of quasar surveys to m, < 22 
will move further down the QLF and increase the density of 
quasars observed, meaning a smaller survey can be used to 
constrain the cluster ing to compar a ble acc uracy as the SDSS 
or 2dF. Using the iHopkins et all d2007bl) QLF to estimate 
the relevant space density of quasars above the flux limit as 
a function of redshift and assuming the errors in clustering 
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FIG. 14. — Characteristic inferred quasar-hosting halo mass corresponding 
to the model clustering as a function of redshift shown in Figure rnn for our 
adopted cosmology. Dotted, d ashe d, and dot-dashed lines are for the appro- 
priate flux limits as in Figure fO] The black, blue, and red lines show the 
left, center, and right models ("efficient feedback", grow th t o typical z ~ 2 
quasars, and "maximal" growth, respectively) from Figure[T5]with effectively 
infinitely deep flux limits (;' < 30); all are identical below z ~ 2.6. 

sume that quasars at a given redshift z > 2, with a typical L* 
and corresponding M BH at that redshift, will grow by the same 
proportional amount as a halo which represents a 6a fluctua- 
tion (characteris tic of halos hosting observed z ~ 6 quasars, 
I Y. Li et ail 120061) from the observed redshift to z = 2, after 
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amplitude relative to those in ICroom et al.l d2005l) scale as 

—l li 

N qS o , we estimate that for the redshift interval 3.5 < z < 4.5 
(3.75 < z < 4.25) a field size of - 25 deg 2 (50deg 2 ) would be 
sufficient to distinguish between the first two models (efficient 
high-redshift feedback and all high-redshift quasars growing 
to z ~ 2 luminosities) at ~ 2.5-3 a. The last model ("max- 
imal evolution") predicts an even more extreme departure in 
clustering properties, and could be distinguished or ruled out 
at ~ 2.5-3 a by clustering observations from 2.75 < z < 3.25 
in just a ~ 8- 15 deg field. 

7. DISCUSSION/CONCLUSIONS 

We compare the clustering of quasars and different galaxy 
populations as a function of morphology, mass, luminosity, 
and redshift, and demonstrate that these comparisons can be 
used to robustly rule out several classes of models for quasar 
triggering and the association between quasar and galaxy 
growth. In each case, the observations favor a model which 
associates quasars with the "formation event" of ellipticals, a 
strong prediction of theoretical models which argue that ma- 
jor, gas-rich merge rs form ellipticals and trigger quasar activ- 
ity (IHopkins et al]l2006bl) . 

The predicted bias as a function of mass/luminosity for sys- 
tems which once hosted quasars agrees well at all masses and 
luminosities with that observed for early-type populations. 
In other words, the clustering of an M g . d \ elliptical galaxy is 
exactly what we would expect if thes e galaxies, which typi- 
cally contain M B h = /xM ga i (m ~ 0.0 01 Mago rrian et al.ll 19981: 
iMarconi & Huntl2003t IHaring & Rixl2004 BHs. represented 
the dominant "hosts" of the quasar population for a brief pe- 
riod, setting at that redshift with an Eddington-limited L* = 
iEdd(^BH) "epoch" of activity. In the most basic sense, this is 
a confirmation that ellipticals today were indeed the host pop- 
ulation of high-red shift quasars (see also Porcian i et al.ll2004 
ICroom et al.l 120051) . with the appropriate corr esponding BH 
masses. This should not be surprising, since the lSoltad(fl982h 
argument demonstrates that most BH mass must have been ac- 
cumulated in bright, near-Eddington "quasar" epochs, and the 
tightness of the local BH-host mas s relation (and sim ilar BH 
mass-host property relations, see iNovak et al.l 120061) argues 
that BH growth must be tightly coupled to the host properties. 
However, there are additional non-trivial implications. 

First, this implies that there really is a characteristic host 
and BH mass "active" at a given epoch, traced by the QLF 
L*. This is an important pr ediction of certain theo retical mod- 
els for quasar lightcurves (Hopkins et al. 2005c d), and sup- 
ported by other lines of observation above. Furthermore, this 
implies that the formation "epoch" for BHs of a given mass 
must be relatively short in time, as continually adding BHs of 
a given mass (at lower Eddington ratio or in radiatively inef- 
ficient states) to the population would dilute the agreement in 
FigureQ] Quasars are active in characteristically different par- 
ent halo populations at different redshifts - i.e. most systems 
cannot undergo multiple separate periods of quasar activity, 
at least at z < 2. We find further support for this by consid- 
ering observed quasar clusterin g as a funct i on of luminosity, 
which favors the predictions of iLidz et al.1 ([2006), namely a 
relatively weak trend of bias as a function of luminosity. In 
fact, the combination of quasar clustering measurements as a 
function of luminosity and redshift supports at high signifi- 
cance previ ous suggestions of little or no luminosity depen- 
dence (e.g., lAdelberger & Steidell [2005al; iMvers et al.ll2006l 
l2007al) . and is inconsistent with the predictions of simpli- 



fied "light bulb" or exponential quasar light curve mod els 
(e.g.. lKauffmann & Haehnell2002l: IWyifhe & Loebll2003l) at 

>4-5<j. 

This relates to a subtle but important distinction: this im- 
plies that the halos of the dominant population of M ga \ ellipti- 
cals are the same halos as those which hosted the correspond- 
ing quasar activity. A significant fraction of the M ga i early- 
type population cannot form from later collapsing halos, as 
this both requires the buildup of BHs of the same mass at a 
different time, ruled out by the observations above, and would 
dilute the clustering agreement. 

Second, the clustering of late-type galaxies at a given lu- 
minosity or mass does not agree with the evolved clustering 
of quasars. This argues that it is specifically the progeni- 
tors of early-type galaxies which hosted quasars. Although 
this is not surprising, given that observations find it is specif- 
ically bulg e mass/velocity dispe r sion which correlates with 
BH mass dGebhardt et all l2000t iFerrarese & Merritt 120001: 
iMcLure & Dunlopll2002t IHaring & Rixll2004l) . there still ex- 
ist classes of models which would generically associate BH 
formation with galaxy formation, star formation, or halo viri- 
alization. Our comparison of observations rules out these sce- 
narios. 

We further invert our comparisons to predict quasar cluster- 
ing as a function of redshift, and compare this with the ob- 
served clustering of red and blue galaxies at each redshift. 
Quasars do not trace a uniform/constant population with red- 
shift - i.e. they are not cosmologically long-liv ed, as has been 
noted in many previous clustering studi es (see iPorciani et alJ 
I2004t lMartinll2004 ICroom et alJl2005l Further, quasars do 
not trace the clustering of "established" red or blue galaxy 
populations. This rules out models in which quasars are asso- 
ciated with cyclic (e.g JCiotti & Ostrikerfl2001blBinneyl2004l) 
or radio "heating" modes over a Hubble time in red galaxies, 
as well as (at least the most straightforward implementations 
of) models which generically a ssociate quasars with star for- 
mation (e.g JGranato et aLll2004 or disk instabilities in high- 
redsh ift, gas-rich disk systems (e.g., iKauffmann & Haehneltl 
l2000l) . 

Note, however, that this does not rule out the presence of 
these accretion modes at low luminosities and/or low red- 
shifts. Many of the observations discussed above limit bright, 
high Eddington ratio quasar activity to a single, short-lived 
epoch. Long-lived accretion in a "radio-mode" is believed 
to be associated with particula rly low Eddington ratio activ- 
ity dHoll2002tlWhite et alj|2006f). perhaps an entirely different 
accretion state (Naravan & Yil ll995: Marche sini etai1l2004 
Uesterl l2005h lPellegrinill200H iKoerding et al.ll2006h . and not 
typical L* QLF activity. Models which invoke "radio mode" 
type accretion at low Eddington ratios in quasar "relics" (e.g., 
ICroton et al.ll2006t iDe Lucia et al.l 120061) are therefore com- 
pletely consistent with the clustering arguments herein. Other 
observation s of black hole mass functions and quasar Edding- 
ton ratios (lYu & Tremalnel l2002t IMcLure & Dunlopl 12004 
iKollmeier et alJl2005l) also rule out "cyclic" models, insofar 
as they attempt to explain more luminous ~ L* quasar ac- 
tivity. However, theoretical models of stochastic, feedback- 
regul ated accretion in gas-rich systems ( Hop kins & Her nquist 
2006) predict that these fueling modes dominate at low red- 
shift and at typical " Seyfert" luminosities at higher redshifts, 
even where mergers (IHopkins et alJl20l^l2007al) may dom- 
inate the bright ~ L* quasar population. We demonstrate 
that the luminosities at which these "Seyfert" accretion modes 
may dominate are sufficiently low that they have no effect on 
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our results. However, we make predictions for future mea- 
surements of clustering as a function of luminosity at mod- 
erate redshifts, which may be able to detect such changes 
in the characteristic host population and fueling mechanisms 
through the cross-correlation of galaxies and faint, X-ray se- 
lected AGN in deep fields at z ~ 1 . 

It is also important to distinguish the processes which 
may be associated with the initial formation of ellipticals 
from their subsequent evolution. Once morphologically trans- 
formed by a gas-rich merger, for example, mass can be moved 
"up" the red sequence (galaxies increased in mass) by gas- 
poor mergers, which will involve neither star formation nor 
quasar activity, but it cannot be added to the red sequence in 
this manner. As noted above, low-luminosity AGN or "radio 
mode" activity, or halo shock "quenching" may be of critical 
importance to suppressing cooling flows and further accretion 
in massive ellipticals. However, these do not appear to be as- 
sociated with the initial formation of an elliptical or triggering 
of traditional, bright, high-redshift quasars. 

It is interesting that, at all redshifts, quasars cluster- 
ing is observed to be intermediate between blue and red 
galaxy clustering. To the extent that halos grow mono- 
tonically with time, intermediate clustering may imply in- 
termediate halo mass and therefore, perhaps, an inter- 
mediate evolutionary stage - i.e. quasars are represen- 
tative of an evolutionary state "between" blue and red 
galaxies. Th is is predicted in many theoretical models 
dGranato et al.ll2004t IScannapieco & Ohll2004t ISpringel et al l 
2005at iMonaco & FontanoJ 120051: IHopkins et al. " 2006b, 
2007al) . which posit that quasar feedback, or a merger- 
triggered quasar phase, is associated with (regardless of 
whether or not it directly causes) the formation of ellipti- 
cals and "transition" of galaxies to the red sequence. Quasars 
do appear to cluster similarly to likely "merger" populations 
(close pairs, post-starburst/E+A galaxies, and sub-millimeter 
galaxies) at both low (z < 0.3) and high (z ~ 2 - 3) red- 
shifts, but more detailed measurements of these populations 
are needed. We confirm the finding from previous studies 
dPorciani et al.ll2004t ICroom et alj|2005t IPorciani & Norbergl 
l2006t Myers et al.l l2006, 2007a) that quasar clustering is con- 
sistent with a constant halo mass ~4x 10 12 h~ l M Q ; inter- 
estingly, similar to the characteristic mass of small galaxy 
groups (at least at low redshift) in which galaxy-galaxy major 
mergers are expected to be most efficient. Further theoreti- 
cal investigation of this, and the possibility that it may man- 
ifest in an excess of quasar clustering on small scales (e.g., 
iHennawi et al.l l2006. but see also Myers et al. 2007b), are im- 
portant subjects for future work. We calculate the implied 
time delay between the star-forming/LBG "phase" of evolu- 
tion, quasar phase, and red galaxy phase, from the evolved 
clustering of halos of a given initial mass. Although this is just 
a "toy" model of self-similar halo evolution, suggestively, the 
time from quasar epoch to "red galaxy" epoch corresponds 
reasonably at all redshifts with the delay expected for stel- 
lar populations to redden to typical red galaxy colors after 
the termination or "shutdown" of star formation. We caution, 
however, that the systematic uncertainties in measurements 
of these clustering strengths (especially at high redshifts) re- 
main a concern, and future studies which compare uniformly 
selected galaxy and quasar populations across a wide range 
of redshifts adopting consistent measurement methods are 
needed to make these conclusions robust. 

We explore this further by considering the age of BHs (i.e. 
the time since the mean "formation" or quasar epoch for BHs 



of a given relic mass) as a function of mass, and comparing 
this with the age of their hosts. The mean stellar age (and 
dispersion in ages) of early-type BH hosts agrees well at all 
masses, implying that quasars are associated with the forma- 
tion epoch of early-type galaxies. Specifically, these are light- 
weighted or single-burst ages of red galaxies, which tend to 
reflect the last significant epoch of star formation - i.e. quasars 
are associated with the last significant epoch (or potentially 
the termination) of star formation in elliptical hosts. This as- 
sociation can be used to accurately predict quasar clustering 
as a function of redshift in a purely empirical manner, with- 
out any assumptions about quasar light curves, lifetimes, or 
Eddington ratios. 

A similar association does not hold for disk-dominated sys- 
tems hosting BHs; i.e. again demonstrating that BH growth 
and quasar activity do not generically trace star formation. 
Likewise, such an association does not hold for dark mat- 
ter halos, meaning that quasars do not generically trace halo 
formation/asse mbly (even accounti ng for the halo downsizing 
effects seen in Neis tein et al.l l2006). They also do not gener- 
ically trace the cross ing of host halos of the c r itical shock 
"quen c hing" mass in iBirnboim & Dekell (120031) : iKeres et ail 
(2005); iDekel & Birnboiml (120061) . Although our compar- 
isons are consistent with the possibility that this halo mass 
does "quench" gas accretion (and thus quasar activity will not 
generally occur at higher masses), there must be a mechanism 
which can trigger quasars before they cross this threshold. 

We emphasize that this does not imply that most of the stars 
in spheroids form in such a short-lived burst contemporaneous 
with their quasar epoch. Direct calculation of the inferred stel- 
lar population ages from line index and SED fitting (following 
Trage r etal]|2000l) for realist ic star formation histor ies from 
the semi-analytic models of lSomerv ille et al. (200Th and hy- 
drodynamical merger simulations of iRobertson et al.l ([2006b) 
suggests that the ages inferred for present early-type galax- 
ies indeed reflect the epoch of the termination of star for- 
mation, even when > 95% of stars are formed over a much 
longer timescale at significantly earlier times (in these cases, 
in qui escent star formation in disks). Indeed, most models 
(e.g. IHopkins et al.ll2006etlCroton et al.l2006tlDe Lucia et all 
2006) predict that only a small fraction of stars in ellipticals 
were formed in merger-induced starbursts. This, combined 
with the lack of a general correlation between star forma- 
tion in disks and quasar activity, supports the hypothesis that 
quasars trace the end of star formation in present spheroids, 
as predicted by models which associate quasar activity with 
mergers (or other mechanisms) that rapidly exhaust gas and 
transform (assemble) disks into spheroids. 

Finally, we extend these empirical predictions of quasar 
clustering to high redshift, and show that the z > 2-3 cluster- 
ing of quasars is dramatically different depending on whether 
or not feedback is efficient at high redshift (i.e. whether or not 
z > 2 quasars "shut down" after their quasar episode, as ob- 
servations show they do at z < 2). Present observations can- 
not distinguish these possibilities, but future quasar samples 
at z ~ 3 -4 with flux limits m, < 22 should be able to break 
the degeneracies. 

Although there are non-negligible uncertainties in these 
conversions to characteristic host masses or luminosities, 
a number of different lines of evidence support their ro- 
bustness. The conversion to a characteristic BH mass 
(and the lack of evolution to z = 0) is determined by 
direct Eddington ratio observations dHeckman et al.l l2004t 
IVestergaardll2004t iMcLure & DunloplfcoO^ iKollmeier et all 
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120051; iHopkins. Hernquist. & Naravanl 120051) and observed 
cosm i c backgroun d constraints (lElvis et aUl2002l: lUeda et all 
2003; Cao ! l2005h . If we instead adopt the simplest em- 
pirically inferred models of quasar lightcur ves from match- 
ing the QLF and local BH mass function (Yu & Tr emainj 
2002tlMarconi et al.l2004l:IShankar et al.l2004tlHopkins et all 
2007b), we come to an identical conclusion. Likewise for 
more sophisticated models which incorporate effects of dif- 
ferent accretion state s and compare radio, X-ray, and optical 
QLFs (Merloni 2004), and theoretical models of quasar light 
curves from numerical hydrodynamic simulations of galaxy 
merg ers which dynamicall y incorporate BH growth and feed- 
back (IHopkins et al.ll2006bl) . Whether we adopt the observed 
BH-host mass/luminosity relations, or those from simulations 
dDi Matteo et al.ll2005HRobertson etal]|2006at IHopkins et all 
l2005allbTlcT) . we obtain an identical result. Adopting just the 
BH-host mass relation, and using the time since the quasar 
epoch to determine the observed M/L (i.e. assuming this rep- 
re sents a stellar age an d using the population synthesis models 
of iBruzual & Charloj (l2003h to predict M/L in a given band) 
also does not change this comparison. 

Finally, direct construction of halo occupation (HOD) 
models from observed quasar clustering dPorciani et all2004i 
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Porciani & Norberd 120061: lAdelberger & Steidell l2005bt 
Croometal 2005)" leads to the same conclusions regarding 
the luminosity L» reflecting the evolving characteristic mass 
of active quasars, their weak subsequent BH growth, and 
host properties at z = 0. In other words, our uncertainties in 
this approach are most likely dominated by the substantial 
measurement errors in the bias b(z) of quasar and galaxy 
populations, not by the systematics in our methodology. As 
such, future improved measurements of quasar clustering and 
bias at high redshifts, particularly as a function of luminosity 
(e.g. using the proximity effect; Faucher-Giguere et al. 
2006), as well as improved galaxy clustering measurements 
which can resolve the clustering as a function of mass or 
luminosity at z > 1, will strengthen the constraints herein 
and continue to inform models for quasar fueling and their 
associations with spheroid formation. 
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